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Abstract — In practical systems, due to the time-varying radio 
channel, the channel state information (CSI) may not be known 
well at both transmitters and receivers. For most of the current 
multiuser multiple-input multiple-output (MIMO) schemes, they 
suffer a significant degression on the performance due to the 
mismatch between the true and estimated CSI. To alleviate 
the performance penalty, a robust downlink multiuser MIMO 
scheme is proposed in this paper by exploiting the channel 
mean and antenna correlation. These channel statistics are more 
stable than the imperfect CSI estimation in the time-varying 
radio channel, and they are used, in the proposed scheme, to 
minimize the total mean squared error under the sum power 
constraint. Simulation results demonstrate that the proposed 
scheme effectively mitigates the performance loss due to the CSI 
mismatch. 

Index Terms — multiuser MIMO, downlink, robust, imperfect 
CSI. 

I. Introduction 

The multiple-input multiple-output (MIMO) system, em- 
ploying multiple transmit and receive antennas, has been rec- 
ognized as an effective way to improve the spectral efficiency 
of the radio channel |[T] ||2| . More recently, multiuser schemes 
have been investigated for MIMO systems to further improve 
the multiuser sum capacity. 

Early studies have assumed a perfectly knowledge of the 
channel state information (CSI) available at the transmitter. |3| 
extended the single-user scheme [4] to the multiuser system. 
However, without exploring the multiuser channel information, 
it simply treated the multiuser interference as the white noise. 
The scheme in |5|, on the contrary, utilized the multiuser 
information effectively to minimize the total mean squared 
error (TMMSE) and, naturally, possessed a better performance. 

The CSI can be obtained at the transmitter either by using 
a feedback channel from the receiver to the transmitter in 
frequency division duplex (FDD) systems, or by invoking the 
channel reciprocity in time division duplex (TDD) systems. 
However, using feedback in FDD systems, the limited re- 
sources for the feedback, associated with the propagation delay 
and schedule lag, heavily degrade the accuracy of the CSI at 
the transmitter As to the channel reciprocity in TDD systems, 
the accuracy of the CSI is corrupted by antenna calibration 
errors and turn-around time delay. In respect that the perfor- 
mance would degrade significantly under the imperfect CSI, 
it is necessary to design a multiuser scheme which is stable 
to the imperfect CSI. 

In robust design methodologies, Maxmin (worst-case) and 
Bayesian (stochastic) are two well known ones Q. The 



former optimizes the performance under the worst case of 
random channels, thus, it is so conservative that its average 
performance is even worse than non-robust schemes |8|. The 
latter maximizes the ensemble average performance over a 
pre-described stochastic distribution of the CSI. When the 
stochastic distribution matches well with the true CSI, the 
latter outperforms the former 

The scheme in [7| was a Bayesian design for downlink 
multiuser MIMO systems with the imperfectly known CSI. 
It introduced a channel error matrix to the cost function of 
101 , then found the solution which minimized the average 
cost. However, similar with f3l, the multiuser interference was 
also treated as the white noise. Therefore, it is expected that 
the performance can be improved by exploring the multiuser 
information. 

In this paper, a robust scheme for downlink multiuser 
MIMO systems is proposed based on the TMMSE criterion. A 
more general channel model involving the channel mean and 
antenna correlation is considered. The scheme is a Bayesian 
design which minimizes the average cost function under the 
sum power constraint. 

The rest of this paper is organized as follows. The channel 
model and problem formulation are described in Section II. 
The Section III presents the design of the robust multiuser 
scheme for the correlated imperfect known channel under 
the sum power constraint. Simulation results and analysis are 
given in Section IV. Finally, the Section V concludes the paper 

Notation: Boldface upper-case letters denote matrices, and 
boldface lower-case letters denote column vectors. tr{-), (•)*, 
(•)^, 1 1 ■ 1 12 and 1 1 • I |_F denote trace, conjugate, conjugate trans- 
position, Euclidian norm and Frobenius norm, respectively. 
E{-) represents the expectation of a stochastic process. 

denote the (i,j)-th element and j-th column of a matrix, 
respectively. 

II. Problem statement 

A. Channel Model 

Consider a base station (BS) with M antennas and K 
mobile stations (MS's) each having Ni{i ^ 1 . . . K) antennas. 
Represented by a matrix Hj e QNixM ^ ^j^g downlink MIMO 
channel to MSi is assumed to be frequency-flat and quasi- 
static block fading. Suppose a non-zero-mean channel with 
both transmit and receive antenna correlations, is written 



as 
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where Wi is the ratio of the power in the mean compo- 
nent to the average power in the variant component of H^; 
A,j G C^i^^i is random, we assume that its entries form 
an independent identical distribution (i.i.d.) complex Gaussian 
collection with zero-mean and identity covariance, i.e., ^ 
CA/'(0, 1); Ho,» G C^^^M jjjg normahzed channel mean, 
and Ro,r.j G C^-^^^ and Rt G C*^^*^ are the normahzed 
correlation matrices of the receiver of MS, and the transmitter 
of BS, respectively. (1) is rewritten into the following for 
simplicity lilli 



(2) 



where = \JWil{Wi + l)Ho,i is the channel mean, and 
R-r,z = + 1)1^0, r.i is the equivalent correlation matrix 

of the receiver of MS^. 

The channel mean and correlation are more stable than 
the instantaneous channel information, and they are usually 
acquired by time-averaging on channel measurements. In the 
Rayleigh channel, for example, the non-zero channel mean 
is obtained by averaging channel measurements over a window 
of tens of the channel coherence time |T0). Furthermore, the 
channel model (2) can also denote the correlated Rician MIMO 
channel, in which case the channel mean represents the line- 
of-sight (LOS) component of the MIMO channel. 

In this paper, we assume that transmitters and receivers only 
know channel means and antenna correlations. 

B. Problem Formulation 

We assume that there are Li{i — \...K) substreams 
between BS and MSi(i = 1 . . . K), that is to say, BS transmits 
Li symbols to MS^ simultaneously. Then the signal received 
at MSi is 

K 

y, - Af H, ^ BfcXfe + Af n, (3) 

fc=i 

where y,; G C^*^^ is the received signal vector, and G 
QLiXi jg jjjg transmitted signal vector from BS to MS^ with 
zero-mean and normalized covariance matrix I. We assume the 
transmitted signal vectors of different users are uncorrected. 



i.e., E (x,x 



= 6ijl, where 6ij is the Kronecker function. 



5ij — 1, when i — j and Sij = 0, when i ^ j. We also assume 
the noise vector is independent of any signal vector. A linear 
post-filter A^ G C^'^^'(i = 1 . . . K) is used at MSi to 
recover an estimation of the transmitted signal vector x^. 
defined in (1) [or (2)] denotes the MIMO channel from BS 
to MS,;. B, G C^^'^^'ii = 1...K) is used at BS to weight 
the transmitted signal vector x,;. After passing through B^, x^ 
becomes into an M x 1 signal vector which is transmitted 
by M transmit antennas of BS. G C^'^'^ is the noise 
vector with the correlation matrix R„. = a'^l^., where In; 
denotes the x N; identity matrix. In this paper, we assume 
Li = ■ ■ ■ = Lk = L. 



A, and B^ (fc = 1 . . . K) are jointly designed to minimize 
the total MSB under the sum power constraint. Hence, we get 



min TMSE^e(y: ||xfc-yfc||' 
Vfe=i 

s.t. tr (^E^BfcBf^ <F 

where P is the total transmit power of BS. 

III. Robust TMMSE Scheme 
According to (3), the j-lh user's MSB is 



(4) 



MSE, 



EiWii, 



K 



= E{tr{Afn,{Y: B,Bf )Hf A, + a^Af A, 
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-BfHfA,-AfH,B,+I)) 

. (5) 

Substitute (2) into (5) and note E (Hi) — H^, hence 
E {iiI a^'R]^ = 0, we obtain 

MSE, = tr (Af H, (g B,Bf )Hf A, + Af A, 
-Bf Hf Aj - Af HjBj + I) + E{tr{Af 
rJ,, A,R? (E B,Bf )Rp Af Rjf A,)) 



(6) 



Observe the last part in (6) 



K 



i?(ir(AfR*,A,R|(EB,Bf)RpAfRjfA,)) 



1 K 



= HRfiE B.Bf )Rp£;(Af R|;A,Af R^A,)) 

2—1 

(7) 

Moreover, as A^ CJ\f{0, 1) with i.i.d. entries, 

E{[A,].^^^{[Aj]. J") = Sn,„a. Therefore, the (TO,n)-th 
entry of the expectation in the right side of (7) is 



iff 



Ei 



AfRlfA.AfRlA, 



m.n 



= tr(Rif A,Af Ri^£;([A,]^ „ ([A,]^,„)^)) 



(8) 



Thus, the expectation in the right side of (7) is 

E{AfRlf AjAfRl^j) ^ tr{Af-RrjA,)I (9) 
Substitute (9) into (7), we obtain 

K 

E 

1=1 

= ir(Rf (£ B,Bf )Rptr(Af R,-,,A,)I) 

i=i 



£;(tr(AfRj^.A,R*(£B,Bf)RpAfRjfA,)) 
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(10) 



Substitute (10) into (6) 



K 



MSEj = tr(AfHj(^ B,Bf)Hf Aj- +a2Af A, 



'■j 



-BfHfA,-AfH,B,+I) 
-ir(Af R,,,A,)tr((E B.Bf )Rt) 



(11) 



The Lagrangian of (4) is 

L(Ai, . . . , Ak, Bi, . . . , Bif) 

= E MSEk + \{tr{ f: BfeBf ) - P) 

k=l k=l 

(12) 

where A is the Lagrangian muhipHer associated with the 
total power constraint. So the Karush-Kuhn-Tucker (KKT) 
conditions [10] of (4) are 



dL{Ai, Ak,Bi, ...,Bk) 
OA* 

dL{Ai,...,AK,Bi,...,BK) 
dB* 



(13) 
(14) 



x(^triJ2B,Bf)-P^ 



^ (15) 
A > (16) 



Among them, (13)(14) come from the fact that the gradients 
of the Lagrangian (12) definitely vanish at the optimal point, 
and (15) is known as the complementary slackness. According 
to (11)--(14), we obtain 

A, = (H,( £ BfcBf )Hf + tr{{ £ B^Bf )Rt)R.,i 

k=l k=l 



+a2l)-iH,B, 



(17) 



B,: = ( E Hf AfcAf Hfc + tr{ £ A^Af R,,fc)Rt 

k=l k=l 

+AI)-iHf A, 

(18) 

Substitute (18) into (15) , we can find A is the root of the 
equation 



A(tr(X(X + Y + AI) ')-P) = 



where 



K 



X = ^HfA,AfHfc 



k=l 
K 



Y = tr{Y^ AfcAf R,,fc)Rt 



k=l 



As Rt is the normalized correlation matrix of the transmitter, 
it is Hermitian, hence both X and Y are Hermitian, and so is 
X + Y. Perform the eigenvalue decomposition 



X + Y = UDU 



H 



where U is unitary and D is diagonal. If A 7^ 0, (19) can be 
rewritten to 

[u^xu]^^ 



^ (dn + A)2 

n— 1 ^ ^ 



p = 



(23) 



where dn is the rt-th diagonal element of D. Using a binary 
search, the root of (23) can be found quickly. Since the left- 
hand side of (23) is monotonous in A when A > 0, the upper 
and lower bounds on A can be acquired by replacing rf„ with 
dmin and drnax, respectively. Thus, 



/ir(X) 



P 



A,, 



/tr(X) 



P 



- d„ 



(24) 



(25) 



where (•)+ means that the expression takes the value inside 
the parentheses if the value is positive, otherwise it takes zero. 
A numerical binary search, then, can be carried out between 
these two bounds to find the root of (23) up to a desired 
precision. Once there is no root between the bounds, which 
implies that the inequality constraint (16) is inactive, A = 
is the only available solution to (18). From (17) (18), it can 
be found that the optimal transmit matrices Bk{k = 1 . . . K) 
are functions of the receive matrices Ak{k = 1 . . .K), and 
vice versa. Therefore an iterative algorithm to calculate Afc 
and Bfc(fc = 1 . . . K) is proposed as follows. 



Initialize ' and A^ ' (k = 1 . . . K) randomly, 
n = 

1) Calculate A from A^"' {k = I . . . K) by 
solving (19). 

2) Calculate B^"+^' (fc = 1 . . . if) from A^"^ 
{k = l...K) and A using (18). 

3) Calculate A^^+^\k = I . . . K) from B<"+'^ 
{k = l...K) using (17) . 

4) Repeat 1), 2) and 3) until 

K 



E(l|Ai"+^'-Ai"'||| 



IB 



(i+l) _ g(")l|2 



fc ll-F, 



< e. 



In our simulation, we set e = 0.0001 . 



(19) 



IV. Simulation results 

(20) In this section, numerical simulations have been carried 
out to evaluate the performance of the proposed scheme. We 
assume that the BS equipped with four antennas {M — 4) is 
communicating with two MS's (K = 2) each with N receive 

(21) antennas (A^i = N2 = N). Also we assume that the number of 
substreams of each MS is equal to 2 (Li = L2 — 2), moreover, 
both the two MS's have the same Wi {Wi ^ W2 = W). 
QPSK is employed in the simulations and no channel coding 
is considered. Let the transmit antenna correlation matrix Rt 
be [Rt] j j — 0.91*^^1 and the receive antennas be uncorrected, 

(22) i.e., R,._i = Iat. 
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Fig. 1 . Comparison of the BER perfomiance of the robust scheme and the 
TMMSE, when N = 2 and W = 10, 50, 200, 1000. 
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Fig. 2. Comparison of the BER perfomiance of the robust scheme and the 
TMMSE, when Af = 2, 3, 4 and = 50. 



Firstly, we compare the bit error rate (BER) of the proposed 
robust TMMSE scheme with that of the traditional TMMSE 
scheme. Defining the signal-to-noise ratio (SNR) as the ratio 
of total transmitted power to the noise power of each antenna 
(SNR = P/afj), Fig. [T] is the average BER curves versus 
the SNR When iV = 2. In order to highlight its impact on 
the BER performance, different values of W are used in the 
evaluations. When W is small, the channel mean poorly re- 
flects the instantaneous channel state, thus the receiver can not 
completely eliminate the interference among the transmitted 
signals, which further induces an irreducible error floor at high 
SNR region. However, the proposed robust scheme overcomes 
the traditional one with a noticeable gain. As the W increases, 
the transmitter obtains more precise CSI, therefore the residual 
interference is mitigated greatly and the error floor vanishes. 
In addition, the gain between the proposed robust scheme and 
the traditional one turns small when the uncertainty of the 
channel state is decreasing. 

In Fig. 121 we compare the BER performance when the 
number of receive antennas is increasing. The additional 
receive antennas provide more spatial diversity gain. In this 
figure, W is fixed to be 50 and N changes from 2 to 4. 
Although both the two schemes explore the additional receive 
diversity gain, the proposed robust scheme obviously has a 
better performance for all A^'s due to its insensitivity to the 
imperfect CSI. 

Fig. [3] shows the average MSE as a function of W, when 
SNR = 20dB and A^ = 2. The average MSE of the proposed 
robust scheme is less than that of the traditional TMMSE 
scheme over all W's. Moreover, compared to the traditional 
TMMSE, the descending slope of the proposed robust scheme 
is flat, which further indicates that its performance is insensi- 
tive to the channel uncertainty. Especially when W becomes 
larger, more reliable CSI is available, therefore closer the two 
curves get. 

Fig. |4] illustrates the convergence property of the proposed 
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Fig. 3. Comparison of the average MSE of the robust scheme and the 
TMMSE, when Af = 2 and SNR = 20dB. 













































































\ 


























































SNR = 


IB 


































SNR = 10 dB 




























SNR = 20 dB 









































2 4 6 8 10 12 14 

number of iterations 



Fig. 4. Convergence of the robust scheme, when N = 2 and W = 100. 



robust scheme, when N = 2, W = 100. The curves of 
the average MSE versus the number of iterations needed 
in different SNR's are plotted. The higher the SNR is, the 
more iterations the proposed scheme runs for to converge. 
Fortunately, for the most SNR's, four iterations are big enough 
to guarantee the convergence. 

V. Conclusion 

In this paper, we investigate a robust Unear processing 
scheme for the downHnk multiuser MIMO system under the 
consideration of imperfect CSI. As the traditional downlink 
multiuser MIMO systems depend on the instantaneous CSI too 
much, they suffer poor performance once the CSI is not accu- 
rate enough. In order to deliver a better performance under the 
imperfect CSI, an iterative Bayesian algorithm which explores 
channel statistics to offer a much more stable description to the 
channel state is developed by minimize the total MSE under 
the sum power constraint. Numerical simulations exhibit the 
proposed robust scheme experiences an obvious performance 
gain over the traditional schemes. In addition, the proposed 
iterative algorithm has a good convergence property - after 
no more than four times of iterations, the algorithm achieves 
convergence. 
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